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To clarify the superconducting gap structure of the spin-triplet superconductor S^RuCU, the in- 
plane thermal conductivity has been measured as a function of relative orientations of the thermal 
flow, the crystal axes, and a magnetic field rotating within the 2D RuC>2 planes. The in-plane 
variation of the thermal conductivity is incompatible with any model with line nodes vertical to the 
2D planes and indicates the existence of horizontal nodes. These results place strong constraints on 
models that attempt to explain the mechanism of the triplet superconductivity. 
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Ever since its discovery in 1994 0, the superconduct- 
ing properties of the layered ruthenatc Sr2Ru04 has been 
attracting a considerable interest. A remarkable feature 
which characterizes this system is the spin-triplet pair- 
ing state with d-vector perpendicular to the conducting 
plane, which has been confirmed by 17 NMR Knight 
shift measurements Moreover, /zSR experiments sug- 
gest that the time reversal symmetry is broken in the 
superconducting state ||. Up to now, the spin triplet 
pairing state is identified only in superfluid 3 He, heavy 
fermion UPt 3 [|, and organic (TMTSF) 2 PF 6 §, though 
it most probably is also realized in the recently discov- 
ered UGe2 [[§. While in 3 He the simplest p-wave pairing 
state is realized, UPt3 seems to be in a more compli- 
cated /-wave state. At an early stage, the gap symmetry 
of Sr2RuC>4 was discussed in analogue with 3 He in which 
Cooper pairs are formed by the ferromagnetic spin fluc- 
tuation. Then the pairing state with the isotropic gap in 
the plane, d(k) — Aoz(k x + ik y ), where Ao is a constant, 
has been proposed as being likely to be realized . 

However, recent experiments have revealed that the 
situation is not so simple. Neutron inelastic scattering 
experiments have shown the existence of strong incom- 
mensurate antiferromagnetic correlations and no sizable 
ferromagnetic spin fluctuation This implies that the 
origin of the triplet pairing is not a simple ferromagnetic 
interaction. Furthermore, the specific heat C p and NMR 
relaxation rate T-j - on very high quality compounds ex- 
hibit the power law dependence of C p oc T 2 |10| and 



T-f cc T 3 jy] at low temperatures, indicating the pres- 
ence of nodal lines in the superconducting gap. These 
results have motivated theorists to propose new models 
which might explain consistently the spin-triplet super- 
conductivity in the ruthenates . Most of them pre- 
dict the line nodes which are vertical to the 2D planes. 
However, the detailed structure of the gap function, es- 
pecially the direction of the nodes, is an unresolved issue. 
Since the superconducting gap function is closely related 
to the pairing interaction, its clarification is crucial for 



understanding the pairing mechanism. 

A powerful tool for probing the anisotropic gap struc- 
ture is the thermal conductivity k, in which only the 
unpaired electrons are responsible for the thermal trans- 
port in the superconducting state. Compared to the spe- 
cific heat and NMR measurements, an important advan- 
tage of the thermal conductivity is that it is a directional 
probe, sensitive to the orientation relative to the ther- 
mal flow, the magnetic field, and nodal directions of the 
order parameter ]Tf^l9|. In fact, a clear 4- fold modula- 
tion of k with an in-plane magnetic field which reflects 
the angular position of nodes of d x i_ y i symmetry has 
been observed in YBa2Cu307_5, demonstrating that the 
thermal conductivity can be a relevant probe of the su- 
perconducting gap structure [ pp|j2l| . Although previous 
attempts have been made to measure the thermal con- 
ductivity in Sr2Ru04, the experimental resolution were 
not good enough to identify the nodal directions . In 
the work reported in this Letter, we have performed a 
high-precision measurement of the in-plane thermal con- 
ductivity as a function of angle between the thermal cur- 
rent q and the magnetic field H rotating within the R11O2 
plane, which is sufficient to resolve the gap structure. 

Several single crystals with different T c 's were grown 
by the floating-zone method. The thermal conductivity 
was measured by a steady state method with one heater 
and two ruthenium-oxide thermometers. In the present 
measurements, it is very important to rotate H within 
the RuC>2 planes with high accuracy because a slight 
field-misalignment produces a large effect on k due to 
the large anisotropy. For this purpose, we constructed a 
system with two superconducting magnets generating H 
in two mutually orthogonal directions and a 3 He cryostat 
equipped on a mechanical rotating stage with a minimum 
step of 1/500 degree at the top of the dewar. Computer- 
controlling two magnets and rotating stage, we were able 
to rotate H continuously within the Ru02 planes with 
a misalignment less than 0.015 degree from the plane, 
which we confirmed by the simultaneous measurement of 
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the resistivity. 

The inset of Fig. 1(a) shows the T-dependence of k/T 
in zero field. Since the electrical resistivity is very small 
which is an order of 0.1/iil-cm, the electron contribution 
well dominates over the phonon contribution |p3| . At the 
superconducting transition, k/T shows a kink. At low 
temperatures, k/T decrease almost linearly with decreas- 
ing T with finite residual values at T = 0. The residual 
k decreases with increasing T c and is very small in the 
crystal with highest T c (=1.45 K). These Independence 
and the residual k/T are consistent with the presence of 
the line nodes 

Figures 1(a) and (b) show the ii-dependence of k for 
the sample with T c =1.45 K in perpendicular (H_L ab- 
plane) and parallel fields (H|| afr-plane), respectively. In 
both orientations, k increases with H after the initial 
decrease at low fields. The consequent minimum is much 
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FIG. 1. Field dependence of the thermal conductivity of 
Sr 2 Ru0 4 (T c =1.45 K) in (a) perpendicular field H± ab 
and (b) parallel field H|| [110]. The thermal current q is 
applied along the [110]-direction. In perpendicular field, 
k is ^-independent below the lower critical field. Inset: 
T-dependence of k/T in zero field for two crystals with dif- 
ferent T c (To = 1.45 K and 1.32 K). 



less pronounced at lower temperatures. At low T, k in- 
creases linearly with H . We note that the f/-linear de- 
pendence of k is observed only in the very clean crystals 
with T c > 1.3 K and k increases with an upward curva- 
ture in samples with lower T c . In parallel field k rises 
very rapidly as H approaches H C 2 and attains its normal 
value with a large slope (d,K/dH), while k in perpendicu- 
lar field remains linear in H up to H C 2- The understand- 
ing of the heat transport for superconductors with nodes 
have largely progressed during past few years |l7],|24|,^5| . 
There, in contrast to the classical superconductors, the 
heat transport is dominated by contributions from delo- 
calized quasiparticle states rather than the bound state 
associated with vortex cores. The most remarkable ef- 
fect on the thermal transport is the Doppler shift of the 
quasiparticle energy spectrum {sip) — > e(p)—v s -p) in the 
circulating supercurrent flow v s P6[ . This effect becomes 
important at such positions where the local energy gap 
becomes smaller than the Doppler shift term (A < v s -p), 
which can be realized in the case of superconductors with 
nodes. In the presence of line nodes where the density 
of states (DOS) of electrons N(e) has a linear energy de- 
pendence (N(e) oc e), N(H) increases in proportion to 
VH. While the Doppler shift enhances the DOS (§6), 
it also leads to a suppression of both the impurity scat- 
tering time and Andreev scattering time off the vortices 
pl| , p5| . This suppression can exceed the parallel rise in 
N(e) at high temperature and low field, which results in 
the nonmonotonic field dependence of k(H). 

It has been shown that in the superconductors with 
line nodes, k increases in proportion to H in the "super- 
clean regime" where the condition, ■£ <C jf^ is satisfied. 
Here Y is the pair breaking parameter estimated from the 
Abrikosov-Gorkov equation *(l/2 + r/27rT c )-*(l/2) = 
ln(T c o/T c ), where ^ is a digamma function and T c o is the 
transition temperature in the absence of the pair break- 
ing. Assuming T c0 =1.50 K and A = 1.76T C , Y/A is es- 
timated to be 0.025 (0.067) for T c =1.45 K(T C =1.37 K), 
showing that our field range is well inside the superclean 
regime except at very low fields smaller than 400 Oe 
(1000 Oe). Thus the iJ-linear dependence of k(H) ob- 
served in very clean crystals is consistent with the k of 
superconductors with line nodes. The steep increase of k 
in the vicinity of H C 2 in parallel field is also observed in 
pure Nb J27j. When the vortices are close enough near 
H C 2, tunneling of the quasiparticle excitations from core 
to core becomes possible, which leads to large enhance- 
ment of quasiparticle mean free path and k. The absence 
of a steep increase in perpendicular field may be related 
to the difference of the vortex core structure. We note 
that a similar behavior is observed in UPt3, in which the 
steep increase of k is present in H|| c while is absent in 
H||6@. 

We now move on to the angular variation of the ther- 
mal conductivity in parallel field. Figures 2(a) and (b) 
depict k(H,0) as a function of 9 = (q,H). No hystere- 
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sis of k related to the pinning of the vortices was ob- 
served in rotating 8. In all data k(H,8) can be de- 
composed into three terms with different symmetries; 
k(H,0) = kq(H) + H2e{H) + m$(H), where k is 6- 
independent, K2b(H) = C^eiH) cos 20 is a term with 2- 
fold symmetry, and K4$(H) = C±e (H ) cos AO with 4-fold 
symmetry with respect to the in-plane rotation. Figures 
3 (a)-(d) show K4s/k u after the subtraction of Ko- and 
K2S-term from k. 

The sign and magnitude of C29 and C40 provide im- 
portant information on the gap structure. The term K29 
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FIG. 2. (a) Angular variation (6 =(q,H)) of K(6)/K n for 
Sr 2 Ru0 4 (T c =1.45 K). q is applied to [110]-direction. (b) 
Same data for the sample with T c =1.37 K. q is applied to 
[100]-direction. The solid lines show the 2-fold component in 
k(6)/k h - 
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FIG. 3. (a)-(d) The 4-fold symmetry K,4g/n n at several 
fields. 



appears as a result of difference of the effective DOS for 
quasiparticles travelling parallel to the vortex and for 
quasiparticles moving in the perpendicular direction. In 
the presence of vertical nodes, the term K^g appears as a 
result of two effect. The first one is the DOS oscillation 
associated with the rotating H within the plane. This 
effect arises from the fact that DOS depends sensitively 
on the angle between H and the direction of nodes of 
order parameter, because the quasiparticles contribute 
to the DOS when their Doppler-shifted energy exceeds 
the local energy gap. In this attains the max- 

imum value when H is directed to the antinodal direc- 
tions and becomes minimum when H is directed along 
the nodal directions Q , E^ . The second one is the the 
quasiparticle lifetime from the Andreev scattering off the 
vortex lattice, which has the same symmetry as the gap 
function |2li| J?l| J25f . This effect is important at very low 
fields where k decreases with H . In addition to the 4- 
fold symmetry associated with vertical nodes, there is 
another contribution to K4g-term, which originates from 
the tetragonal band structure inherent to the Sr2Ru04 
crystal. We will discuss this effect later. 

The most important subject is "Is the observed K40 a 
consequence of the vertical line nodes?" Before analyzing 
the data, we list up the various proposed gap functions 



1. Type-I: Vertical nodes at 
|pHl5|; d(k) = A z(s'mk x - 

A[)zk x ky(k x + iky}. 



(±tt,0) and (0,±tt) 
-ismky) and d(k) = 



2. Type-II: Vertical nodes at (±tt, ±tt) jj||j2§]; d(k) = 

A z{kl-kl){k x +ik v ). 

3. Type-Ill: Horizontal nodes 0,^; d(k) = 
Aoz(k x + ik y )(cosck z + a) with a < 1 (c is the 
intcrlayer distance) and d(k) = Aozk z (k x + ik y ) 2 . 



As shown in Figs.3(a)-(c), K40 shows minimum at H | 
[110]. Therefore, this result immediately excludes the type 
I symmetry, in which k±q should exhibit a maximum at 
H || [110]. We next discuss the amplitude of K40. Fig- 
ure 4 depicts the fZ-dependence of |C 2 g| and \C^g\. In 
the vicinity of H c2 where K increases steeply, \C4g\/K n 
is of the order of a several % (see Fig. 3(a)). However, 
|C4g|/K n decreases rapidly and is about 0.2-0.3% at lower 
field where k increases linearly with H (see Figs. 3(b) 
and (c)). At very low field where K decreases with H, no 
discernible 4-fold oscillation is observed within the reso- 
lution of \C ie \/K n < 0.1% (see Fig. 3(d)). 

Recently, the amplitudes of /C40 for various symme- 
tries with vertical nodes have been calculated at the 
field range where k obeys an H- linear dependence [15|. 
We will examine our result in accordance with Ref. [15|. 
For both type I and II symmetries, \Cie\/K n is expected 
to be about 6% at low field. Apparently, the observed 
\C±e\/K n < 0.3 % at low fields are less than 1/20 of the 
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prediction for type I and II symmetries. Thus it is very 
unlikely that the observed J^-fold symmetry is an indica- 
tion of vertical line nodes. We then consider the tetrag- 
onal band structure as an origin of n^g. This effect can 
be roughly estimated by the in-plane anisotropy of H c2 
pOf . In our crystal, we find that H c2 is well expressed as 
H c2 ((t>)/H c2 (0) = 1 + A cos 40 with A = -0.013, where 4> 
is the angle between H and a-axis. In Fig. 4, we plot 
| C40 1 = \A\Hdn(H)/dH calculated from the in-plane 
anisotropy of H c2 with no fitting parameter. The cal- 
culation reproduces the data, indicating that the 4-fold 
symmetry of n is indeed mainly due to tetragonal band 
structure. 
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We next discuss K 2 g which provides an additional im- 
portant information on the gap structure. According to 
Ref. H, a large 2-fold amplitude, \C 28 \/K n > 25% is 
expected for type I and II symmetries when q is in- 
jected parallel to the nodes. To check this, we applied 
q along [110] and [100] directions as shown in Figs. 2 (a) 
and (b). In both cases \C 2 g\j K n is about 1 %, which is 
again much less than expected for the case of vertical line 
nodes. Thus both 2- and 4-fold symmetries of the ther- 
mal conductivity are incompatible with any model with 
vertical line nodes. 

We now examine the Type-Ill symmetry without n^g- 
term associated with the nodes. The magnitude of C 2 $ 
provides a clue toward distinguishing between the two 
gap functions listed under the category of type-Ill. Ac- 
cording to Ref. fid), a large magnitude of \C 2 g\/K n > 



30% is expected for d{k) — Aozk z (k x + iky) 2 ■ In fact, 
a large 2-fold oscillation is observed in the B-phase of 
UPt3 with this symmetry |l6|. On the other hand, 
much a smaller \C 2 q\/ K n ~8% is expected at T = 
for d{k) = Aoz(k x + ik y )(cos ck z + a). Although the 
value is still several times larger (which may be due to 
finite temperature effect which reduce |C 2 © | ) , it is much 
closer to the experimental result. These results lead us 
to conclude that the gap symmetry which is most con- 
sistent with the in-plane variation of thermal conduc- 
tivity is d(k) = Aoz(k x + ik y )(cosck z + a), in which 
the substantial portion of the Cooper pairs occurs be- 
tween the neighboring Ru02 planes. These results im- 
pose strong constraints on models that attempt to ex- 
plain the mechanism of the triplet superconductivity. We 
finally comment on the orbital-dependent superconduc- 
tivity scenario, in which three different bands have dif- 
ferent superconducting gaps ||. In this case, our main 
conclusion can be applicable to the band with the largest 
gap (presumably the 7-band). 

In summary, the in-plane thermal conductivity of 
Sr2Ru04 have been measured in H rotating within the 
planes. The angular dependence is incompatible with any 
model with vertical line nodes and strongly indicated the 
presence of horizontal line nodes. 
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